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LOW-LEVEL WATER VAPOR HELDS FROM THE 
VISSR ATMOSPHERIC SOUNDER (VAS) *miT W1 
CHANNELS AT 11 AND 12 MICRONS 


D. Ch«sten, L. Uccellini, tad W. Robinson 


ABSTRACT 

A stiiss of higli-resolution water vapor fields are derived fliom the 
1 and 1 2 ttm channels of the VISSR Atmospherk Sounder (VAS) on 
GOES-S. The low4evel tropo^heiic moisture content is separated 
from the surfKe and atmo^hezk radiances by using tlw differaitial 
absorption across the **split window** along with the average air temper- 
ature frmn imbedded radiosondes. Fields of predpitable water (gm cm*^) 
are presented in a time sequence of five false color images taken over 
the United States at 3-hour intervab on 13 July 1981. Vivid subsynop- 
tic and mesoscale patterns evolve at 1 S km horizontai resolution over 
the 12-hour observing period. Convective cloud bimations develop 
from several areas of enhanced low-level water vapor, especially where 
the vertical water vapor gradient (inferred from overlays of high-level 
moisture patterns found in the corresponding VAS 6.7 fim channel) is 
relathrely strong. Independent verification at radiosonde sites indicates 
fairly good absolute accuracy (±0.67 gm cm’^ from 1.8 to 5.3 gm cm*^), 
and the spatial and temporal continuity of the water vapor features 
indicates very good relative accuracy. Residual erron are dominated by 
radiometer none and unresolved clouds. The algorithm can produce 
quality controlled images of the vertically integrated water vapor in 
relatively cloud free areas from real-time VAS data, closely monitoring 
the development of regional moisture features which would remain 
otherwise undetected and/or unresolved. 
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LOW-UVEL WATER VAKMl HELDS FROM THB 
VISSR ATMOSWIERIC SOUNDER (VAS) *miT WINDOW” 
CHANNELS AT 11 AND 12 MICRONS 


1.0 INTRODUCTION 

1.1 Mothntion 

Th« VISSR Atmospheric Sound«Hr (VAS) is the thennal infrutd ndicmieter now opertting on 
the (jeostationaiy Operational Environmental Satellite (GOES). The VAS instrument is an un- 
proved version of the previous imager on riie GOES, with eleven infrared chanmls added to the 
1 1 ttm window and visible detectors. Radiometric quality has been ui^ded as well, wirii ten bit 
precision in place of eight, ±1.5”K calibration accuracy instead of t3.0*K, and six infrared detectors 
instead of one. The VAS is intended to exploit its geosynchronous station by observing the devel- 
opment of mesoscale temperature and water vapor fields with frequent multiH:hannel images at 
1 S km horizontal resolution (nadir view) in the spectral windows and molecular absorption bands 
of the earth’s thermal infrared spectrum (SucMni, et al., 1971). In particular, the pair of channels 
at 1 1 and 12 Mm are designed to use the differential water vapor absorption across this part of the 
spectrum in order to estimate the water vapor absorption along the line of right. The two channels 
are called a ”split window” because they use the difference in brightness at two adjacent bandpasses 
in order to separate the atmospheric and surface radiance contributions. The VAS split window is 
used in a simple algorithm to produce false color images of the water vapor content of the lower 
troposphere. This bulk water vapor estimate is certainly not as sophisticated as a water vapor pro- 
fUi retrieval using all of the VAS channels and human quality control. Rather, it is a real-time 
algorithm for making high-resolution mesoscale images from only two channels which display the 
horizontal low-level water vapor distribution in relatively clear air. 

For the near future, NOAA will operate VAS with minimal impu:t upon the pre-established 
community of VISSR (visible and 1 1 Mm infrared) data users. The additional VAS detectors make 
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it possible to deliver the 6.7 and 1 2 iim radiances simultaneously along with the conventional 
VISSR data to any ground station with VAS data handling capability. Consequently, these three 
VAS infrared channels can be used for routine mesoscale monitoring of water vapor fields. These 
expectations have motivated the development of an algorithm which can automatically produce 
real-time water vapor estimates from these channels alone. 

1.2 Background 

In general, atmo^heiic water vapor sounding from satellite radiances is more difficult than 
temperature sounding because; 

a. Vertical water vapor structures are often not resolved by the broad vertical extent of the 
passive radiometer response functions. 

b. Most of the precipitable water along a line of sight is concentrated in the air just above a 
bright background of highly variable topography, emissivity and skin temperature. 

c. The radiance contribution from the atmo^here is a combination of both the air tempera- 
ture and the water vapor content 

Nevertheless, some quantitative infonnation about the atmospheric water vapor can be retrieved 
from passive thermal infrared satellite data. 

For example, recent case studies of mesoscale water vapor retrievals with the High-Resolution 
Infrared Sounder (HIRS) over the Texas-Oklahoma area were able to display dryline features down 
to the 30 km horizontal resolution of that instrument (Hillger and Yonder Haar, 1981). These 
HIRS analyses could only determine the relative horizontal structure of the water vapor features, 
and not their absolute magnitudes. Similar error patterns have also been found in the operational 
satellite temperature soundings: there are larger errors near the surface, biases in areas of unusual 

topography or clouds, and systematically underestimated horizontal gradients. Such problems 
have been found both over ocean (Phillips, et al., 1979) and overland (Schlatter, 1981). Likewise, 
simulated case studies of VAS water vapor retrievals within pre-convective environments have 
shown a pattern of good relative gradient determination, but poor absolute accuracy (Chesters, et ai.. 



1982). In particular, the VAS simulations showed a poor linear correlation between the split win- 
dow brightness difference and the low-level mixing ratio, due to the genwal problems with water 
vapor retrievab. 

The first VAS instrument was launched on GOES-4 in the autumn of 1980, and it is currently 
operational at 13S‘W. Initial experience with the VAS data indicated that the instrument is a rela- 
tively stable, somewhat biased radiometer (Menzel, et al„ 1981). Temperature and water vapor 
profiles have been retrieved by using careful radiometric corrections and a physical retrieval method 
at hand-picked, clear flelds-of-view (Smith, et al., 1981). The VAS on GOES-4 has also been used 
as a mesoscale temperature and water vapor field interpolator between radiosonde stations (Chesters, 
et al., 1981). Subsynoptic low-level water vapor features definitely appear in objective analyses of 
the VAS soundings from both retrieval methods. Mesoscale features derived from VAS cannot be 
verified with the relatively coarse, infrequent radiosonde network. In this sense, VAS is the only 
available source of frequent, quantitative, high-resolution atmospheric data for the United States. 
The data used in this paper are from the second VAS instrument, launched on C}OES-S in the spring 
of 1 98 1 . The data were acquired during demonstration and testing of the VAS Processor at NASA/ 
GSFC in the summer of 1981, Just before GOES-S was moved from its post-launch station at 8S^W 
to its current operational station at 75^ W. 

1 3 Outline 

This paper describes a fast, simple algorithm for estimating the low-level precipitable water in 
clear air from the VAS split window channels at 1 1 and \2fim. The low-level moisture algorithm 
presented in Section 2 is based on a physical model for radiation transfer through a single thick 
layer of moist air in the lower troposphere. The two-channel algorithm involves tliree unknowns: 
surface temperature, air temperature, and integrated water vapor. Consequently, the calculation 
requires one empirical parameter derived from conventional weather data. We use the “average 
brightness temperature" for the lower troposphere, determined at clear radiosonde sites imbedded 
within the VAS radiance field. The algorithm is refined on simulated VAS data in order to 



partmetexize the tvenge molecular abaoiptioR and to establiah quality controls. Section 3 demmi* 
stntes the low-level moisture algorithm with one day of VAS data. A false cokir image sequence 
of precipitable water fields over the United States is computed from the split window radiances. 

The low'^evel water vapor patterns are also compared to the upper-4evel water vapor patterns infer- 
red from the VAS 6.7 pm radiances. Areas are noted in which VAS indkates a strong vertical water 
vapor gradient before the outbreak of convective instabilities. The error budget for the split win- 
dow alg<mthm is analysed, and the sateUite-derived precipitable water esthnates ate verified by 
comparison to independent radiosonde reports. Section 4 summarizes the resulte and discusses 
possible extensions of the low-level moisture a^orithm. 

2.0 MODEL FOR LOW-LEVEL MOISTURE IN THE VAS SPLIT WINDOW 

2.1 The VAS Split Window , 

A “split window** requires two channels at adjacent wavelengths in a fairiy transparent region 
of the earth’s atmospheric spectrum. Because the radiances in both channels are subject to nearly 
the same modifications by surface emissivity, aerosols and clouds, the difference in brightness tem- 
perature is controlled by the differential molecular absorption. The molecular and surface contribu- 
tions to the radiance can be separated algebraically, and the molecular concentration can be 
determined from calculated values of tlie absorption cross sections. Naturally, the accuracy of this 
method depends upon several factors, including the accuracy of the radiometer, the field-^f-view 
registration between the channels, the parameterization of molecular absorption cross sections, and 
the formulation of an algorithm which is robust enough to operate with only two channels and 
simple quality controls. 

The VAS split window covers the 1 1 to 13 Mm region of the earth’s radiance spectrum, sketched 
in Figure 1. The window region is bracketed by the 10 ^m ozone band and the 15 Mm carbon 
dioxide band, with increasing molecular absorption at the longer wavelengths due to water vapor 
and carbon dioxide. The differential absorption was originally proposed as the basis for a correction 
to the sea surface temperature derived from the VISSR 1 1 Mm window (Prahabrkara, et al.. 1974). 
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By contrast, our algorithm removts the background skin temperature in order to cakulate the 
water vapor absorption in the lower troposphere. 

The standard radiance weighting functions dr/dln? are shown in Figure 2 for the VAS 11,12 
and 6.7 channels, and their shapes reflect the vertical water vapor profile of the U.S. Standard 
Atmo^here. The 1 1 and 1 2 im cluumela transfer radiance from the entire lower tropoq>here (600 
to 1 000 mb), with the radiance exchange concentrated in the boundary layer, where most of the 
water content occurs. The other moisture-sensitive channel is centered on the strong 6.7 ttm water 
vapor absorption band indicated in Figure 1 . This bandpass has already been tested from polar 
orbit with NIMBUS-4 data (Rodgen, et al, 1976) and from ^synchronous orbit with METEOSAT 
data (Ramond, et al., 1981). The high-level 6.7 ftm water vapor patterns can be compared to the 
low-level 1 1 and 1 2 tm water vapor patterns in order to infer the vertical water vapor structure 
(Petersen, et al., 1982). 

2.2 The Single-Layer Model 

A single-layer radiation model for the low-level water vapor is suited to the limited vertical 
resolution of the VAS 1 1 and 1 2 Mm channels shown in Figure 2. The main features of the model 
are: a surface radiating at an effective temperature T^, a layer of air radiating at an average temper- 
ature T,, and a satellite-observed brightness temperature of T*,, in channel v. The water vapor- 
determined transmissivity controls the blend of surface and air brightness in a linearized radiation 
transfer equation: 


T*y * Tj Ty + Tg ( 1 - r^), 

- (T% - Tg)/(Tg - Tg). 


( 1 ) 


In the fust form of ( I), surface radiance is attenuated by the factor r and is replaced by the 
corresponding fraction 1 - r from the air. The net effect is controlled by both the surface-to-air 
temperature contrast and the available water vapor. For example, Tg * 280*’K air attenuating the 
view of a Tj * 300*K surface produces a T* * 292*K scene when r ■ 0.60. Previous VAS studies 
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(Chesters, et al., 198 1 end 1982) attempted to use the simple algebraic difTerence between the ^lit 
window brightness temperatures, T* ^ * measure of the water vapor content. However, 

both theory and practice show that a considerable fraction of the underiying surface temperature 
confuses this water vapor indicator. 

In the second form of ( 1), the transmissivity is expressed as the Action of the surface-to-air 
temperature contrast which is visible from the satellite. The transmissivity ratio between the 1 1 
and 12 Mm split window then removes the surface temperature algebraically: 

n:"ll * - w„ - T,). (2) 

The steps required for the application of (2) to the VAS split window are; (a) a simple parameteri- 
zation linking the transmissivity ratio to the water vapor content, and (b) an empirical method for 
determining the average air temperature. 

2.3 Parameterizing the Transmissivity 

Transmissivity liv, P) is the probability' of the satellite observing a photon at wavenumber v 
after it has passed up through the atmosphere from pressure level P. The total transmissivity can be 
computed exactly from radiosonde temperature T(P) and water vapor Q(P) profiles as the product 
of the dry transmissivity through the well mixed gases (CO 2 , NO, CH^, N 2 and O 2 ) times the wet 
transmissivity over the water vapor pro^e: 

^tot , ^dry^wet 

t(v, P) ■ exp(-k(*', P)sec(d)), (3) 

/* Spice 

Q(P')da(»;.P',T(P')). 

Exact transmissivity calculations like (3) must be done laboriously from molecular line-by-line 
cross sections (McClatchey, et al., 1973) and p-type and e-type continuum cross sections (Roberts, 
et al., 1976). whose sum provides the absorption cross section da (v. P'. T(P')) at pressure P'. The 
absorption coefficient kti' P> is the vertical integral over the product of the temperature-sensitive 
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cxx>n section da (v, P\ T(P')) and the molecular concentration profile Q<P')< The absorption along 
a line of sight is proportional to the projection factor sec(9), where 9 is the zenith ai^ from a 
fleld-of-view to the satellite through a plane parallel atmoq)here. Minor effects due to ozone ab- 
sorption, the carbon dioxide continuum, and aerosol scattering are nei^ected in our calculations. 

For the VAS bands, the water vapor absorption cross section is approximately twice as great 
at 1 2 Min as at 1 1 Mm. The transmissivity in the VAS 1 1 Mm channel is typically r|j ■ 0.60. Con- 
sequently, the satellite receives roughly 60% of the photons emitted by the surface at 1 1 Mm, and 
only 36% (Le., 0.60^) at 1 2 Mm. For practical use in real-time computatioiu, the transminrvity at 
every VAS field-of-view is estimated from average values of the wet and dry absorption for the 
single layer atmoq>heric model. 

The water vapor absorption at 1 1 and 1 2 Mm is parameterized as the product of an average 
absorption cross section a^^^ and the precipitable water PW (measured in gm cm*^) projected on 
the iine>of-sight'. 

^wet M exp(--k'^**sec(d)), 

fcwet.awetp^ (4) 

PW - g-1 / ^ CXP) dP. 

•'surface 

The absorption coefficient for the well mixed gases is parameterized with a linear term to 
account for the modest e^ect of air temperahire upon the CO-, absorption: 

jdry . exp sect®)), 

^ ( 5 ) 

kdry . k + (dk/dT) (T^ - 280). 

The average wet and dry absorption parameters (a^^^ k and dk/dT) in (4) and (5) are determined 
by least squares fits to exact multi-layer transnissivity calculations for a w-ie ranging (arctic to 
tropical) set of 32 radiosonde reports, simulated at zenith angles from O'* to 60”. The average values 
are listed in Table 1 for the VAS 1 i and 12 Mm channels. They should apply with acceptable 


Table 1. 

Wet and dry absorption parameters for the lower tix>posphere, deiiml fnMn least square fits 
to wide ranging simulations of trananisshrity for the q)lit window channels at 1 1 and 1 2 Mm. 


Channel 

V (cra*^ ) 

anh 

k 

dk/dT(K"l) 

12 Mm 

789.24 

0.3169 

0.06114 

0.00091 

11 Mm 

897.40 

0.1591 

0.01066 

0.00019 


precision to any clear VAS field-of-view. The RMS error over the global range of simulated condi- 
tions was ±0.035 in the transmissivity ratio ^ 

VAS data is approximately ±0.025, corresponding to only ±0.3 gm cm'^ uncertainty in a moisture 
burden of 5.0 gm cm~^ along the line-of-sight. Consequently, errors due to the transmis^ty 
parameterization will contribute less than 1091) uncertainty to the satellite-derived water burden. 

2.4 The Average Air Brightness Temperature 

Finally, precipitable water is related to the brightness temperatures and average air temperature 
by combining (2), (4) and (5) and the information in Table 1 into the following set of equations: 


wet / _ wet 
12 /^l 


(T*12 - - V 

expi-Ak'^ sec(ff)) 


( 6 ) 


Ak^^y - (*tl2-*^ll> ^ (dkjj/dT-dkjj/dT)(Tj-280K). 

Given the satellite observed bilghtnesses, T*j2 and T*j j . one can first use f6) to estimate from 
PW values at some radiosonde sites imbedded within the radiance field. Then, one can use this 
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avenge ndiosonde>based value of to make satellite estimates of PW om the rest of ^e fMd. 

The model represented by (6) actually treats the “average air tempenture** as an empirical 
parameter which accounts for the lower tropospheric radiance. An even simple empirical method 
for determining Tg is to adjust it by hand until the satellite-derived precipitable water values nearly 
match the values expected within the scene. This human interactive method is especially useful 
when the scene’s precipitable water content can only be estimated from climatologkal con»lma- 
tions, such as for oceanic fields. Of course, more elabonte horizontal models for Tg could be 
devised, such as gradients, air masses or objective analyses between sites. 

The sensitivity of satellite-derived precipitable water to the air temperature can be estimated 
from (6) as the derivative dPW/dTg. For typical mid-summer conditions of T^ « 300 K, Tg ■ 280 K, 
PW » 4 gm cm“^ and sec ($) * 1 25, the senativity is dPW/dTg * 0.6 gm cm*^ K“^ . The senativity 
to radiometric errors is a comparable value. Such sensitivity is undoubtedly relate to the substan- 
tial absolute errors which were found in previous studies of satellite-derived water vapor fields. 

Like those studies, this report finds evidence for very good determination of water vapor gradients, 
but less certain absolute accuracy. Indeer* the average air temperature Tg in (6) acts like a water 
vapor scaling parameter, and the empirical approach to its determination at ground truth sites is 
intended to cancel as many biases as possible. 

3.0 VAS SPLIT WINDOW OBSERVATIONS OF LOW-LEVEL WATER VAPOR 

3.1 Case Study, 13 July 1981 

At 1200 GMT on 13 July 1981, a weak surface low was drifting slowly from eastern Colorado 
toward Nebraska, with a cold front extending Lorn Nebraska toward northern Illinois and into the 
southe.Tt Great Lakes region (Figure 3-A). The southeastern United States was dominated by a 
broad high pressure system centered near the Mississippi-Alabama border. The resultant anticyclo- 
nic circulation forced a southerly air flow over a rather broad band west of the Mississippi River. 

This air was noticeably and uniformly humid with surface dewpoints greater than 20*C over the 
region from the Gulf of Mexico northward into southern Minnesota, and with dewpoints greater 
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than 24^C located both in southern Illinois and near the Gulf Coast Precipitation observed by 
radar (Figure 3*^) was mainly confined to thunderstorms along and to the north cf the cold front, 
extending along a line from South Dakota to the Ohio valley. Weak thunderstorms ai'.d showers 
were also reported in the Texas-Loutsiana coastal area. 

By 18(X) GMT, the southerly air flow and the morning insolation had forced surface tempera- 
tures to rise well above 30"C over the central United States, whfle the dewpoints remained high 
(Figure 3-B). Two distinct dewpoint maxima are observed, one directed from the Gulf of Mexico 
toward eastern Oklahoma and the other extending from Iowa to Ohio. A noticeable gradient in the 
dewpoints is evident in the area between eastern Colorado and eastern Kansas. The thunderstorms 
originally in South Dakota had drifted northward and diminished in intensity, while the storms 
originally in Iowa had weakened considerably as they moved into Illinois (Figure 3-E). The convec- 
tion along the Gulf Coast continued to intensify, and numerous rain and thundershowers were 
reported east of the Mississippi River. 

By (KXX) GMT on 14 July 1981 , the surface dewpoints exceeded 24^C in two regions: in the 
Gulf States region extending northward to eastern Oklahoma, and immediately ahead of the now 
warm front in Iowa and west central Illinois (Figure 3-C). The dewpoint gradient from eastern 
Nebraska to Kansas strengthened considerably as drier air moved east-northeast from Colorado into 
Nebraska. Scattered rainshowers had developed in the area of Kansas, Oklahoma and Arkansas, and 
heavier thunderstorms developed from New Mexico to western Nebraska, while the baiid of thun- 
derstorms continued along the Texas-Gulf Coast region (Figure 3-F In Iowa, intense thunder- 
storms with tops approaching 1 8,000 m developed south of the frontal zone and were accompanied 
by weaker cells in extreme eastern Nebraska. 

Our application of the VAS split window will focus upon the changes in the water vapor Held 
in the middle portion of the United States, where (a) convection persisted along the Gulf Coast, 

(b) scattered rain showers developed in Oklahoma, and (c) intense thunderstorms went through 
various stages of development in Iowa and Illinois. 
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3.2 Satellite Imafeiy and Quality Control 


Satellite images also help to set the scene for application of the louMevel water v^K>r algo- 
rithm. Figure 4 contains four satellite images of the United States at midHiioming on 13 July 
1981: 

a. Figure 4-A is an operational visible VISSR image taken from SMS-3 at 1430 GMT. At 
1 km resolution, it shows that the central United States is relatively cloud free, demon- 
strating that unresolved clouds are not a serious hazard at this time to proctsuag the VAS 
thennal infrared data at 1 5 km resolution. Visible imagery is not available at 1 200 GMT, 
due to the low sun angle at that time, so that we must assume sunilar cloud conditions 
during the preceeding VAS and radiosonde observations. Viabte imagery taken after 1 500 
GMT shows the development of broken clouds over much of the east-central United Sta^ 
along with the development of overcast conditions in the regions of precipitation found in 
Figures 3-E and 3-F. 

b. Figua 4-B is the VAS 1 1 ^m infrared image taken from GOES-5 at ISOO GMT. For inter- 
channel comparisons, the radiances have been converted to brightness temperature over a 
range from 200” to 320”K, with a false color spectrum assigned in 10”K steps. Each color 
is shaded internally to display detailed gradients more clearly. A typical clear pixel in the 
central United States has T*|j between 290” and 300”K, with false color red. One can see 
a 10”K gradient between the high plains and the central river valleys, but still little contrast 
between the land and water areas. By 2 100 GMT, the VAS 1 1 ^m images of the central 
United States brightened by another 10”K, while the clear oceanic areas remained about 
the same. This is a substantial dynamic range of underlying surface temperatures to remove 
from the split window water vapor estimates. The cloudy areas in Figure 4-B are radio- 
metrically colder, from 220” to 280” K, and are colored with yellows, greens, blues and 
whites. Radic*onde sites which were judged to be cloud free at 1 200 GMT arc indicated 
with their station names on this 1 500 GMT image. 
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c. Figuie 4-C is the corren>on(Ung VAS 12 fun imtge at 1500 GMT. SinM the w^ervapor 
attenuation is greater at 1 2 Min. the brightness tempeiitures are S* to 10*K km than fouiKt 
in the 1 1 Mm window. The central plains show T* J 2 between 280* and 290*K, shifted by 
almost one 10*K false color step (from red to orange) with re^>ect to Figure 4-B, wlule 
the unattenuated cloud top temperatures remain the same. 

d. Figure 4>D is the cotre 4 >onding VAS 6.7 Mm bri^tness temperature image at 1 500 GMT. 
As expected from the radiance westing function in Figure 2, this channel appean much 
colder than the window charmds. The surface features are totally hMden by the upper- 
tropo^heric water eapor. and only the tops of tlw h^er clouds iq[>p^ as ^landmarks** 
which can be found in the other panels of Figure 4. Suteynoptte upper4evel water rapor 
patterns t^pear as cold, relatively wet areas (colored dark green, near 240*K) mixed with 
warm, relatively dry areas (colored light green and yellow, i^ar 250*K). The warmest 
streaks are interpreted as relatively transparent “dry slots** in the upper air, often associated 
with upper>level jets. A more detailed analysis of the upper-level water vapor patterns 
inferred from the VAS 6.7 fua channel has been reported separately for this case (Petersen, 
et al. 1982). 

The S to 10*K difference between T*j| and T *^2 m Figures 4-B and -C is used as the satellite 
“signal** for precipitable water estimates. 

The noise in the VAS split window channels is approximately ±0.5*K (Chesten. et al.. 1981). 
Since the split window brighmess tempera tuie difference is 5* to 10*K. the signal/noise ratio is 
approximately 1 0/ 1 . Consequently, the water vapor estimates are noticeably limited by radiometric 
noise, even in clear air with a strong surface-to-air temperature contrast. When the T* difference is 
small or when either channel approaches T^. the PW estimate from (6) becomes unreliable. Conse- 
quently, a series of simple thte^old tests is used to identify doubtful flelds-of-view: 

Reject ifT* — Tj< l*Kin either chaimel, 

or T*j 2 - T*n < I*K, (7) 

or PW > 10 gm cm"^. 
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The air temperature threshold accounts for most of the rejected pixels, i^kh are the cold clmid 
tops found in Figure 4. The two-channel threshold reacts pixels where the water vapor is covered 
by low clouds or fog. The precipitabte water threshold is set rather high in order to provide some 
subjective quality control. Because PW greater than 6 gm cm*^ is physkaUy unreasonable for our 
case, the excessive PW valun indicate the location of clouds which are unresolved by the VAS 1 5 km 
footprint, warning the user to regard those areas with suspicion. 

# 

3 3 Empirical Air Temperature Estimate 

The 1200 GMT lower tropo^here’s average air brightness temperature T^ was found by using 
radiosonde-derived PW values in (6) and (7) at the eleven clear fields-of-vkw indicated in Figure 
4->B. The radiosonde-based air temperatures had no systematic variations within the 1200 GMT 
field, so that the aven^ value (282.7 ± 1 3^K) was adopted. For this midsummer case, the average 
value of Tj^ was nearly equal to the average radiosonde temperature reported at 700 mb. Because 
characterizes the mean brightness of the lower troposphere, its value should be little affected by 
dhimal heating of the boundary layer. Consequently, the value of T^ from the 1200 GMT radio- 
sondes was extended to the VAS split window fields throughout the day, up to and including the 
subsequent radiosonde launch at 00(X) GMT. 

Of course, a single value for the air temperature cannot apply over an extremely large domain. 
For instance, the mean value of T^ from the central United States will produce both overestimates 
of the precipitabte water in the distinctly colder air mass over the Great Lakes region and under- 
estimates in the subtropics. The simple threshold tests in (7) will reject the substantially colder air 
mass in Canada (note the low 1 1 nm brightness temperatures in Figures 4-B and ~C) and will even 
reject some of the cold waten in the Great Lakes themselves. 

3.4 False Color Images of Low-Level Water Vapor 

Figure 5 is a series of five false color images of the low-level water vapor derived from the VAS 
split window. The algorithm in (6) was applied to every pixel which was not rejected by the quality 
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controls in (7). The images were taken at thxee-hour intervals from 1200 GMT on 13 July 1981 to 
0000 GMT on 14 July 1981. The PW range is fh)m 0 to 10 gm with a change of hue on the 
false color spectrum every 1 gm cm*^ . Each hue is internally shaded to display detailed PW varia- 
tions, Fields-of-view which were rejected by the quality control thresholds are colored black, and 
most of them obviously corre^ond to the cold cloud tops seen in Figure 4. Drier pixels are colored 
red, normal pixels are colored orange and yellow, wetter pixels are colored light green, and exces- 
sively wet pixels are colored dark green and blue. 

For quality assurance, a histogram has been made of the precipitable water distribution of the 
processed pixels within each frame of Figure S. The number of pixels with PW between 0 and 
10 gm cm~^ is counted within 0.1 25 gm cm'^ bins. The pixels with PW greater than 6 gm cm~* are 
probably contaminated by unresolved clouds. While they account for about 20% of the processed 
fields-of-view, they are clearly grouped into cloudiike patterns within the Gulf of Mexico and along 
the edges of the recognizable cloud masses. The histogram statistics remain nearly stationary over 
the 1 2-hour observing period, supporting the assumption that the average 1 200 GMT air brightness 
temperature applies consistently throughout the day. There is a small shift toward lower average 
PW values by midday, which is consistent with a small acUial increase in T^ due to daytime heating 
of the boundary layer at the bonom of the lower troposphere. 

The low-level water vapor patterns in Figure S show vivid subsynoptic and mesoscale features 
at 15 km resolution evolving over the ! 2-hour observing period. The algorithm has eliminated all 
traces of the underlying 20** K range in skin temperatures observed in the 1 1 pm radiances during 
the day. The spatial and temporal continuity of the features is convincing qualitative evidence of 
the relative accuracy of the algorithm. For instance, the relatively cloud free area in the north- 
western Gulf of Mexico serves as a “secondary standard” for the low-level water vapor, monitoring 
the temporal stability of the algorithm within a region of steady atmospheric and surface conditions. 
Likewise, the dry area in Nebraska monitors the temporal stability of the satellite-derived low-level 
water vapor within a region of steady atmospheric conditions but rapidly varying surface 
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temperature. A fUll spectrum of horizontal moisture variations is evident in Figure 5. Qualitatively, 
the amplitude of the variations is rather mild on the synoptic scale (due to the extension of a slngte 
value for the air temperature to the entire domain), quite large on the subsynoptic scale, mo<terate 
on the mesoscale, and again low on the local scale, with very L^^e pixeMo-pixel variation which is 
not due to unresolved clouds. 

On the subsynoptic scale, the VAS water vapor features in F^re 5 corre^nd to the general 
dewpoint patterns in the conventional nuface analyses in Figure 3. For instance, the moist cloudy 
area in Arkansas drifts slowly into Oklahoma with the low-level flow, aiul develops convective 
elements by 2KX> GMT. The surface dewpoint maximum in the Gulf States is somewhat masked by 
clouds in the satellite predpitabls water images. Subsynoptic details include the intensifying arc of 
clouds and water vapor which propagates from the Gulf of Mexico into Texas. Likewise, the mois- 
ture gradient from Nebraska to Kansas (surface dewpoints from 10^ to 20°C in Figure 3-A) is de- 
tected by the VAS split window (precipitable water from 2 to 5 gm cm“^ in Figure 5-A). This 
gradient strengthens in both analyses over the 1 2-hour observing period, as a markedly dry feature 
moves northeastward across Nebraska, growing to the south of a developing low pressure area. In 
fact, all of the VAS-derived subsynoptic moisture features over the central United States are con- 
firmed at the lower resolution of the conventional analyses. The regional moisture features which 
are normally unobserved in the conventionally data-void areas, such as the Gulf of Mexico and the 
Atlantic Ocean, are subjectively verified by the continuity of their development. 

On the mesoscale, the VAS split window images in Figure 5 display vivid structures and signi- 
ficant changes within the subsynoptic features. For instance, the satellite-derived water vapor 
images depict the evolution of the broad moisture maximum in Arkansas and Kansas into two 
separate moisture bands within which convective clouds ultimately develop. By 0000 GMT, the 
precipitation pattern in the Oklahoma region (Figure 3-F) reflects the distinct low-level moisture 
pattern (Figure 5-E). Likewise, Iowa was the scene of interesting changes during the 1 2-hour 
period. After the morning thunderstorms dissipate over eastern Iowa, a water vapor maximum re- 
mains at 2 100 GMT (Figure 5-D). By 0000 GMT. new thunderstorms develop in southern and 
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eastern Iowa, within a local water vapor maximum and along a water vapor gradient, relatively 
(see Figures 3-F and S-E). The storms along the water vapor gradient in southern Iowa are located 
150 km behind the warm frontal zone, which remains in northern Iowa during this period (Figure 
3-C). Another interesting example of mesoscale variability is found within the stable dry feature in 
Nebraska. The single dry zone originally in southwest Nebraska slowly propagates northeastward 
between 1 500 and 1 800 GMT and then is rapidly enlarged by the development of a second dry zone 
in northeastern Nebraska between 2 1(X) and 0(XX) GMT (Figures 5-G, and •£). Thus, the in- 
crease in the conventional surface dewpoint gradient across eutem Nebraska (Figures 3«^ and -C) 
appears to be the remit not only of the advection of the original dry air but also the rapid forma- 
tion of a second dry region, perhaps due to the downward tran^ort of mid-tropospheric air. 

Figure 5 reveals a wealth of high-resolution details and rapid changes within the VAS-derived 
low-level moisture fields which cannot be fully confirmed by the conventional radiosonde or sur- 
face networks. The changing mesoscale moisture fields in the Iowa region provide an excellent ex- 
ample of the limitations of using conventional analyses for verification of satellite data. There are 
no radiosonde stations in the entire state of Iowa to corroborate the VAS-observed precipitable 
water in the lower troposphere, and the hourly surface data lacks the vertical information necessary 
to verify integrated water vapor content. In particular, the conventional surface stations within 
Iowa and central Illinois report dewpoints approaching 24°C in a band extending immediately south 
of the warm front (Figure 3-C), where the satellite measurements of precipitable water are also 
large. However, the surface dewpoints are relatively uniform throughout Iowa, Illinois and Missouri, 
while the VAS precipitable water field (Figure 5-D) exhibits a narrow band of significantly drier air 
which extends from central Nebraska through northern Missouri, and into Illinois. Therefore, the 
VAS split window apparently can differentiate between those i reas in which the water vapor ex- 
tends over a deep layer (more able to supply convective cells) and those areas in which the water 
vapor is confined to a rather shallow boundary layer (less able to support convection). Another 
example of the mesoscale limitations of conventional analyses is revealed by comparing the smooth 
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surface dewpoint gradient extending from Nebra^a to Oklahoma with the complex VAS moisture 
patterns in the sair^ region. The VAS integrated water vapor images show a separate band of mois- 
tuK and showers from north Texas to central Kansas at 0000 GMT (Figure S-E) which is located 
to the west of the dewpoint gradient analyzed from the surface data (Figure 3<-C). Apparently, 
moist air from the east is overrunning a shallow dry layer which still envelops the surface-4evel 
stations. Once again, convective clouds and rain developed within a VASnletermined water vapor 
maximum, as they did in Iowa. 

In summary, the sequence of low>level water vapor fields in Figure 5 demonstrates that the 
VAS instrument meets its objectives of monitoring rapid mesoscale variations from geosynchronous 
station and of filling the space-time gaps in the conventional regional coverage. 

3.5 Contrast to High-Level Water Vapor 

The radiance weighting functions in Figure 2 indicate that the low-level water vapor patterns 
from the 1 1 and 1 2 channels are radiometrically independent of the high-level patterns observed 

in the 6.7 fim channel, so that a correlation between the two images reflects an actual physical cor- 
relation in the vertical water vapor structure. For example, the 1500 GMT images in Figures 5-B 
and 4-D show similar patterns in the moisture structure over the Great Lakes region, but different 
patterns in the moisture structure over the central Great Plains. In fact, the 0000 GMT thunder* 
storms in Iowa were preceeded by the advection of dry air (inferred from the 6.7 fim channel) over 
the low-level moisture concentration (computed from the VAS split window). A detailed compari- 
son of the entire series of low- and high-level water vapor images is the subjec. of a separate paper 
(Petersen, et al., 1982). 

3.6 Quantitative Verification and Error Analysis 

For quantitative verification of the satellite-derived, low-levcl water vapor values, the conven- 
tional radiosonde network is the only available source of “ground truth” data. The 1 200 GMT 
reports on 13 July 1981 were used to establish the average air temperature, so they do not provide a 
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comptetely independent measure of the error. At that time, from the appearance of the VAS im- 
ages, from radiometric quality controls, and from the stations’ o^ n reports of cloud cover, eleven 
sites were judged to be cloud free. The open circles in Figure 6 compare the satellite-derived pre- 
cipitable water to the radiosonde-derived values, showing an RMS difference of ±0.88 gm cm~^. 

The erron are dominated by a disagreement of 1 .94 gm cm*^ at International Falls, Minnesota, 
which is located within a different air mass from most of the other training sites in the central 
United States. 

The 0000 GMT 14 July 1981 radiosonde reports provide a comptetely indeptndent quantita- 
tive measure of the accuracy, twelve hours after the determination of the air temperature. At this 
later time, twenty-four sites were accepted by the algorithm s threshold tests. The filled circles in 
Figure 6 compare the satellite-derived precipitable water to the radiosonde valuer. The correlation 
coefficient of the independent set is 0.78, with an RMS error of ±0.67 gm cm~^ over a range of 1.8 
to 5J gm cm”^. The mean error for the 0000 GMT verification sites was -0.28 gm cm"^, which 
could be explained by an actual inci«ase in the average air temperatur? Hy ■H).4^K during the twelve 
daylight hours following the 1 200 GMT determination of T^. The wont error in the verification is 
associated with Sault Sainte Marie, Michigan, which reported overcast conditions at 0000 GMT and 
which is located in a different air mass from the 1 200 GMT training area in the United States. Some 
of the other test sites were also reporting scattered or broken clouds at 0000 GMT. Obviously, the 
threshold tests on the split window do not remove all cloud contamination from the water vapor 
estimate. 

The quantitative verification of the satellite-derived precipitable water at “ground truth” sites 
indicates fairly good absolute accuracy over a large range of water vapor content. The qualitative 
verification provided by the spatial and temporal continuity of the low-level water vapor fields indi- 
cates very good relative accuracy at high horizontal resolution. The error budget is dominated by 
radiometric noise, clouds and model limitations: 


18 



a. The VAS ndtometric eiron tie approximately t0.5*K compared to a typical 5.0*K con- 
trast in the split window brightness temperatures. The radiometric noise contribute about 
±0.4 gm cm*^ to the error budget for this case. 

b. Unresolved clouds cause systematic errors whkh are more difHcult to quantify. Simula- 
tions with cloud cover terms in ( 1 ) indicate that clouds colder than the average air temper- 
ature cause precipitable water overestimates in (6), and warm clouds or fog cause under- 
estimates. Unfortunately, the use of more severe thresholds for the split window would 
suppress genuine subsynoptic features in the colder areas, such as the Great Lakes region. 

c. The errors due to simplifications by the linearized sin^e layer model for the lower tropo- 
sphere are also difficult to quantify. From simulations, the random errors in the absorption 
coefficient parameterization are estimated to contribute approximately ±0J gm cm*^. 

The systematic errors from using a simple average air temperature are scarcely detectable 
(less than 0 J gm cm~^) in the statistical, horizontal and temporal behavior of the precipi- 
table water fields in Figure 5. Unfortunately, a more ccMnplicated model would introduce 
more unknowns to be determined. The use of more VAS channels would likewise intro- 
duce more unknowns, by sacrificing the split window’s cancellation of several variables in 
the fleld-of-view, such as mrface emissivity. 

Altogether, the random radiometric and parametric erron contribute approximately ±0.5 gm cm'^ 
uncertainty to the precipitable water estimate, accounting for one-half of the variance in the 
observed ±0.7 gm cm*^ RMS enor. The remaining variance is due to more systematic errors from 
unresolved clouds, model limitations and air temperature variations. 

In summary', the quantitative and qualitative verification of low-level water vapor from the 
VAS split window si.ows fairly good absolute and very good relative accuracy, respectively. The 
algorithm is robust enough to extract precipitable water information from just two VAS window 
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channeU. The background surface temperature vaiiationa are eliminated quite effecthrt ly, and the 
residual errors are dominated by noise and unresolved clouds. 

4.0 SUMMARY AND DISCUSSION 

Fields of low-level atmospheric water vapor are derived from the **split window" 1 1 and 
1 2 Mm infrared channels of the VISSR Atro^heric Sounder ( VAS) carried on the GOES-5 satel- 
lite. The differential water vapor absorption between the two channels allows the separation of 
surface and atmospheric radiances within a fleld-of-view. The average air brightness tenrtperature. 
determined at radiosonde sites, allows the further separation of atmospheric emission and absorp- 
tion termi A simply parameterized algorithm is developed from a single-layer model for the split 
window radiances. The vertically integrated water vapor content (measured in gm cm'^) of the 
lower troposphere is then estimated from a time-sequence of five VAS split window images taken 
every 3 hours over the United States on 13 July 1981. False color images of the integrated water 
vapor show vivid subsynoptic and mesoscale details at 1 5 km horizontal resolution which evolve 
continuously during the 1 2-hour period. Precipitable water estimates are verified with independent 
radiosonde measurements taken 1 2 hours after the initial air temperature determination, and an 
accuracy of ±0.67 gm cm"^ is established over a range of 1.8 to 5.3 gm cm"^. Goudy flelds-of- 
view are detected as being colder than the air, and are simply rejected from the images. Unresolved 
clouds are readily identified as cloudlike patterns of excessive water vapor in the images. Several 
areas of enhanced low-level water vapor develop into convective cloud formations. Convective 
cloud formation is most vigorous in an area of initially clear air with a strong vertical water vapor 
gradient (where the upper-level water vapor content is inferred directly from the VAS 6.7 nm radi- 
ances). The sequence of low-level water vapor fields in Figure 5 demonstrates that the VAS instru- 
ment meets its objectives of monitoring rapid mesoscale variations from geor/nchronous station and 
of filling the space-time gaps in the conventional regional coverage. 

This midsummer day of VAS data provides a wealth of detail for satellite retrieval and meso- 
scale development studies. Vertical water vapor structure has already been inferred by comparing 


:o 


the 6.7 Mm chinnel to the low-level moiitiiie fiehlt derived fxom the qdit whidow, end a |»«limiii- 
aiy meteorolt^ical analyris hat been made frmn meaoacale fiekU of VAS mulri-channel amuM^^ 
(Petersen, et al., 1982). The split windo^v estimates are currently being comi^red to water vapw 
profile retrievals derived from all the VAS chaiuiels and other conventional data, and a detailed 
meteorological analysU is being made of the VAS-derived mesoscaie developments associated with 
the Iowa thundersionns. 

The split window, low-level moisture algorithm itself will be tested and extended with more 
case studies. For instance, the two window channels will be supplemented with the VAS low-level 
temperature sounding channels in order to provide a direct estimate for the “air brightness tempera- 
ture'* at each fteld-of-view, instead of using an empirteal avenge from imbedded radiosonde sites. 
Likewise, the other infrared channels can provide better quality control over the cloud-contaminated 
flelds-of-view. In addition, collocated visible data is available at ! km resolution in the VAS data 
stream which could be used to detect clouds that are not resolved with the 15 km thennal detectors. 
The data-void oceans might be processed more accurately with a rearrangement of the algorithm to 
use the sea surface temperature instead of the average air temperature as ancillary conventional 
data. The algorithm will probably be less effective in the drier winter air over the United States, in 
the cloud- and water vapor-obscured tropics, and in any situation where the surface-to-air temper- 
ature contrast is small. Although automatic processing of two VAS channels is not as accurate as 
human interactive, multi-channel retrievals followed by an objective analysis, the real-time image 
sequences can draw a user’s attention to otherwise undetected mesoscaie water vapor features within 
relatively clear air. Ultimately, the utility of the VAS-derived water vapor patterns will be judged 
by their acceptance as a mesoscaie nowcasting and forecasting tool, when the VAS moisture chan- 
nels are distributed as a product of the operational GOES system. 
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FIGURE CAPTIONS 


Figure 1 : The VAS water vapor channels illustrated by (A) a sketch of the earth’s radiance ^MCtnim 
compared to (B) the fQter req>onse functions at 1 1, 12 and 6.7 ^tm. 

Figure 2: Radiance weighting functions for the VAS water vapor channels at 1 1 , 1 2 and 6.7 fan. 

Figure 3: NWS maps of conventional surface data shown at: (A) 1200 GMT on 13 July 1981, (B) 
1800 GMT on 13 July 1981 and (C) 0000 GMT on 14 July 1981; the solid lines are iso- 
bars (e.g. 16 » 1016 mb) and the dashed lines are isodrosotherms (**C). The ccme^nd- 
ing radar precipitation maps are shown at : (D) 1 135 GMT, (E) 1 735 GMT and (F't 
2335 GMT on 1 3 July 1981; maximum tops (dkm) are underlined, and an arrow with a 
number indicates direction and ^eed (m s*^) for the convective system. 

F^re 4: Midmoming satellite images of the United States on 13 July 1981 : (A) the VISSR visible 
image from SMS-3 at 1430 GMT, (B) the VAS 1 1 lua image from GOES-5 at 1 500 GMT 
(color coded as brightness temperature from 200 to 320 K, and marked with radiosonde 
sites which were cloud free at 1 200 GMT), (C) the corresponding VAS 1 2 {im channel 
and (D) the corresponding VAS 6.7 ^m channel. 

Figure 5 : Image sequence of five low-Ievel water vapor estimates from the VAS split window over 
the United States, made at 3-hour intervals starting from 1200 GMT on 13 July 1981. 

The images are encoded with a false color spectrum for precipitable water values from 0 
to 10 gm cm~^. Pixels below the assumed air temperature are assumed to be cold cloud 
tops and are colored black. A histogram of the number of pixels with satellite-derived 
PW values is displayed beneath each image, with a distribution that remains nearly con- 
stant throughout the day. 

Figure 6 : Plot of radiosonde-derived “ground truth” precipitable water values compared with the 
corresponding VAS-derived values: (o) 1 1 clear radiosonde sites at 1 200 GMT on 13 
July 1981 , used to determine T 3 for the scene throughout the day, and (•) 24 independ- 
ent radiosonde sites at 0000 GMT on 1 4 July 1981. 
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